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The local order around alkali (Li + and Na + ) and alkaline-earth (Be + , Mg + and Ca + ) ions in 
4 He clusters has been studied using ground-state path integral Monte Carlo calculations. We ap- 
ply a criterion based on multipole dynamical correlations to discriminate between solid-like versus 
liquid-like behavior of the 4 He shells coating the ions. As it was earlier suggested by experimental 
measurements in bulk 4 He, our findings indicate that Be + produces a solid-like ("snowball") struc- 
ture, similarly to alkali ions and in contrast to the more liquid-like 4 He structure embedding heavier 
alkaline-earth ions. 
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I. INTRODUCTION 

In past decades, the implantation of ions in liquid he- 
lium has often been used to probe the bulk properties 
of the superfluid stated The attention has also focused 
on the details of the structure of the helium solvent 
around the ions^^ As a consequence of electrostriction, 
a strong increase in the helium density with respect to its 
bulk value is expected in the surroundings of the impu- 
rity sitei^ Within a crude picture, ions in liquid helium 
can be divided in two categories, on the basis of the dif- 
ferent local structure they create in the solvent. Strongly 
attractive ions tend to form a solid-like structure, the so- 
called snowball, characterized by a very inhomogeneous 
helium density in the surroundings of the impurity. Al- 
kali ions are believed to belong to this category. In con- 
trast, singly charged alkaline-earth cations, due to their 
larger radii and weaker interaction with the He atoms, 
are expected to produce a cavity surrounded by com- 
pressed, less structured, and most likely not solidified, 
helium {bubble in the following) £ 

Several investigations of the ion-helium structure have 
been based on measurements of ions mobilities in bulk hc- 
lium££ In particular, the experimental results of Ref. 
have shown that the temperature dependence of the mo- 
bility of Bc + is much more similar to that of alkali ions 
(K+,Rb+ and Cs+)£ and He+£ than to that of the other 
alkaline-earth ones (Mg + ,Ca + ,Ba + ,Sr + )j£ It was thus 
suggested that Be + forms a snowball, similarly to He + 
and alkali ions, and at variance with heavier alkaline- 
earth ions, where "bubble" states are believed to develop 
instead^ The recent discovery that positive ions can be 
captured in helium nanodropletsi promises an extension 
to charged particles of experimental techniques which 
proved extraordinarily fruitful for neutral species.— 

Quantum many-body simulation methods allow for a 
microscopic description of such systems ; 11 ^ 2 as well as for 
the investigation of solid-like behavior in quantum clus- 
ters^ Using quantum Monte Carlo (QMC) simulations, 
within the shadow wave function approach, different pos- 



itive ions have been studied both in bulk 4 He , n ' 14 and 
4 He clusters^ with the aim of computing the effective 
masses of charged impurities moving in liquid 4 He and 
of establishing, at the same time, the surrounding local 
structure of the solvent (for a calculation of the effective 
masses within the hydrodynamical model sec Ref. El). 

The distinction between solid- and liquid-like behavior 
of 4 He around atomic/molecular impurities is a subtle 
issue which requires a careful examination of the dynam- 
ical evolution of the system. Within reptation quantum 
Monte Carlo (RQMC)^ii 7 - the use of suitable (imagi- 
nary) time-correlation functions has been exploited to 
propose a dynamical criterion for discriminating between 
freezing versus quantum melting in small para-hydrogen 
clusters.— In the present paper, we resort to the same 
criterion to investigate the structural order in 4 He clus- 
ters doped with alkali and alkaline-earth positive ions, 
simulated by a more standard ground-state path inte- 
gral (PIGS) Monte Carlo algorithm^ We take Li+ and 
Na + as representatives of alkali metals, while among the 
alkaline-earths, we consider Be + , Mg + and Ca + . Lithium 
and sodium -doped 4 He clusters have been chosen be- 
cause the strong interaction of these ions with the He 
solvent makes them prototypes of snowball structures, 
whose expected solid-like signature can provide a useful 
benchmark in the study of the other ions, where a clear- 
cut distinction between bubbles and snowballs might not 
be easy to be resolved. As a paradigmatic case of liq- 
uid helium- impurity structures, we also simulated a small 
4 Hc cluster seeded with a carbon monoxide molecule. We 
considered ion impurities in clusters rather than in bulk 
helium in order to avoid the possibility of finding ficti- 
tious structures imposed by the symmetry of the period- 
ically repeated simulation box. Moreover, the local 4 He 
structure around the ion in large enough clusters should 
be very similar to that developed in bulk 4 He>^ 

We show that the criterion proposed in Ref. [HI is ap- 
plicable also to larger systems than those investigated in 
that paper, and that it is able to discriminate between a 
crystalline or a liquid behavior of the impurity-He struc- 
ture. In qualitative agreement with experimental data£ 



2 



our simulations show in particular that the 4 Hc structure 
around Bc + is much more similar to that found around 
alkali ions, where a snowball structure is found, than to 
the one around heavier alkaline-earth ions, where 4 He 
exhibits a more liquid-like behavior. 

Our paper is organized as follows. In Sec. II we briefly 
review our theoretical framework. Section III A contains 
our results for very small 4 He clusters (seeded with Li + 
and CO), that have been considered for testing our com- 
putational strategy. In Sec. Ill B we present simulations 
for several ions solvated in larger 4 He clusters, and we 
discuss the implications of these results on the main is- 
sue we are addressing in this paper (i.e. the local struc- 
ture of the 4 He solvent around Bc + in comparison with 
the findings for alkali ions and other alkaline-earth ones). 
Finally, Section IV contains our conclusions. 



II. THEORY 

Our clusters are described by a realistic Hamiltonian, 
H , in which the N 4 He atoms and the ionic impurity are 
treated as point particles interacting via pair potentials. 
When the CO molecule is used as dopant, we neglect its 
internal degrees of freedom and we treat it as a rigid lin- 
ear rotor with only translational and rotational degrees 
of freedom. The interactions are derived from accurate 
quantum chemistry calculations, either by an interpola- 
tion of tabulated values, as for Li + -He and Ca + -He ) 19 i 20 
or by parametrized analytical expression, as for He-He* 21 - 
He-CO^ and for the remaining He-ions pairsi 23 i 24 The 
size of the ion-doped clusters varies up to a maximum of 
70 4 Hc atoms, whereas, in the case of CO, we considered 
a smaller cluster, with only fifteen particles (COQHcis), 
approximatively corresponding to one complete solvation 
shcll.-^ 

We simulate the system using the ground-state path in- 
tegral Monte Carlo scheme^ where the imaginary-time 



evolution operator, 



is exploited to improve sys- 



tematically a suitable trial function ^t, by projecting it 

out onto the state \typ) = e -/3ff I^t)- The latter con- 
verges to the exact ground state, ^f , for (3 — > oo. By 
Trotter slicing the imaginary-time evolution into m time 
steps of length e = (3/m, and by choosing an appro- 
priate definite positive short-time approximation for the 

imaginary-time propagator, G(R',R;e) tts (R'|e _eH |R) 
between the configurations R (a vector which describes 
the coordinates of all 4 He atoms in the cluster, in addi- 
tion to those of the impurity) and R', expectation values 
of quantum operators, O, on the state can be ex- 

pressed in a path-integral form: 



(O) = 



(*p\0\*p) 



dXO{K m )P(X), (1) 



system, and O is diagonal in the coordinate represen- 
tation. By sampling the probability P(X) via a gener- 
alized Metropolis algorithm, and accumulating the val- 
ues assumed by the operator O over the sampled paths, 
the desired expectation value can be accessed within a 
known statistical error and without any mixed estimate 
nor population-control bias. 

For the relatively small CO@Hci5, it is sufficient to 
use the primitive approximation^ for the imaginary-time 
propagator G, and then to sample the probability P(X) 
by the RQMC algorithm^ 1 ! In the case of ionic im- 
purities, we resorted to the pair approximation for 
employing the bisection-multilevel sampling procedure.— 
The dependence of (0)p on (3 and e needs to be checked, 
in order to extrapolate the result to the exact value. Our 
projection times j3 varied between 0.1K -1 and 0.7K -1 , 
while we used time steps e ranging between 0.0025 K _1 
and 0.0125 K -1 for clusters seeded with ions, and equal 
to 0.001 K" 1 for CO@Hei 5 . These values of (3 and e have 
been empirically found to ensure the stability of our re- 
sults. 

We used a Jastrow-type trial wave function: 



t = cxp|^ - y^ui(rj) - y^M 2 (^ ) 
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where P(X) = * T (Ro)IlS G(R i ,R i _ 1 ;e)* r (R 2TO ), 
X = {Rq, • • ■ , R2m} is a path of configurations of the 



where rj and indicate, respectively, the distances be- 
tween the z-th 4 He atom and the ion, and that between 
the i-th. and the j-th 4 He atoms. The radial functions 
Hi and U2 have been optimized independently for each 
one of the simulated systems, by minimizing the cluster 
variational energy with respect to a small number (10-12) 
of adjustable parameters. For COQHcis, the one-body 
term U\ also depends on the angle between the molecular 
axis and the vector position of the 4 He atom with respect 
to the molecule center of mass (details can be found in 
Ref. [H). 

Within the path integral scheme, O can be ei- 
ther a static operator (such as the Hamiltonian 
itself) or a time-dependent one, such as those 
needed to calculate imaginary-time correlation functions, 
Acxp(— tH)P>. Their evaluation requires an extended 
time range, 2(3 + t, but it is nevertheless straightfor- 
ward: cab{t) = (^p\A(to)B(to + r)\^p) is estimated 
as Si^( R i)- B ( R i+fc)): where k = r/e, (■) indicates 
a statistical average over the sampled paths, and the n% 
values of i span the range (3/e < i < 2m — 2(3 /e — k. 
The results that will be shown in the following Sections 
have been obtained by varying r within 1 K _1 . Notice 
that time-translation invariance implies that these corre- 
lations are independent of tq . 

As it was proposed in Ref. [l3l . imaginary-time correla- 
tions can be exploited to unveil either a solid- or a liquid- 
like behavior in quantum clusters, by monitoring the per- 
sistence of geometrical signatures of a solid- like structure. 
The criterion of Ref. [ID is based on the definition of the 
rotationally invariant time-correlation functions 
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FIG. 1: Left panel: Radial density distribution of 4 He atoms 
around the Li + ion in Li + @He]v for N = 8 (solid line), N = 9 
(long dashed line), and N = 10 (short dashed line). Right 
panel: 4 He chemical potential pjv = En — En-i as a function 
of the cluster size, N, in Li + @Hejv, whose total energy is En- 

„ 'Em(Qlm( t o)Qlm( t o + r)) , , 

C M T ) = „ , — i—-^ , , (oj 

of the multipole moments, 

Q LA/ = VM(2LTT) / drp(r)r L F L * M (0, 0), (4) 

of the cluster mass distribution p(r) , around its center of 
mass. In Eq. [31 Qlm(tq) is the multipole moment of the 
system at (imaginary) time To and Qlm( t o + t ) indicates 
the multipole calculated after the configuration at time 
To + t has been rotated back in order to minimize the 
particles diffusion between the two time instants. 

For a rigid system cj,(t) is independent of time, 
whereas it decays to zero at large r if the system un- 
dergoes random fluctuations. A large value of cl(t) at 
large times is thus a distinctive mark of the persistence 
of the shape ("rigidity") of the system. This criterion 
has been originally applied to small clusters of para- 
hydrogen molecules (both pure an doped with a CO 
molecule), for cluster sizes around twelve molecules, at 
which a first solvation shell is completed As a pre- 
liminary step in the study of the correlation functions 
Cl[t), we evaluate the static multipole correlation func- 
tion, qL = J2m(Q*lmQlm), which provides a rotation- 
ally invariant characterization of a rigid bodyJ^ This 
quantity is not suitable for discriminating the cases where 
a given multipole vanishes on average but q^ is different 
from zero because of fluctuations, from those where a 
non- vanishing q^ value is due to the average shape of the 
system. However, non-vanishing q^ values allow to se- 
lect a set of L's among which to search for slow-decaying 
time-correlations Cl(t). 
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FIG. 2: Upper panel: Multipole time-correlation functions 
cl(t) of Li + @He8 (L — 4: diamonds, L = 5: filled circles) 
and CO@Hei5 (L = 5: triangles). Lower panel: Diffusion 
coefficient of the 4 He atoms in Li + @Heg, as calculated from 
the original configurations (empty circles) and from configu- 
rations rotated back as described in the text (filled circles). 



III. RESULTS AND DISCUSSION 

A. Testing the approach on small clusters 

As a test case, we first discuss the results for a Li + 
ion in small 4 He clusters (Li + @Hejv). The left panel of 
Fig. [1] displays the radial distribution functions of 4 He 
atoms around the impurity ion, for cluster sizes in the 
range N = 8 — 10. For > 9 the magnitude of the max- 
imum of these functions stays nearly constant while the 
tail of the distribution extends to larger distances from 
the dopant. Correspondingly, the chemical potential rises 
steeply from N = 8 to N = 9. These findings indicate 
that for N = 8 a first solvation shell has been completed. 
We notice that an even more marked jump in the chem- 
ical potential occurs at N = 6, thus suggesting a great 
stability of Li + @He6- In particular some of the mul- 
tipole imaginary-time correlation functions of Li + @Hcg 
show an extremely slow decay in imaginary time, which 
is indicative of the rigidity of the system (see below). 
However, since we aim at revealing the existence of solid- 
like signatures within the first shell of larger clusters, we 
will not discuss any longer the case of Li + @Heg, focus- 
ing instead on the solvent structure in a cluster of eight 
4 He atoms, i.e. one representing a fully developed first 
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FIG. 3: He density distribution for Li+@He 8 . The dis- 
played isosurface contains about 50% of the total number of 
He atoms. The Li + ion (not shown) is located at the center 
of the cluster. 



solvation shell. 

The upper panel of Fig. [2] shows the multipole time- 
correlation functions of Li + @Hes for L = 4, 5 and, for 
comparison, of CO@Hei5 for L = 5. All the multi- 
poles correlations of CO@Hei5 decay rapidly, indicating 
a liquid-like behavior of the first He shell around the 
molecule, and we arbitrarily chose to show the L = 5 
curve as representative of the general behavior of the 
cl(t) correlations for any L. On the other hand, in 
Li + @Heg, all the multipole time-correlations with L < 4 
decay rapidly, whereas for L = 4 and, more evidently, 
for L = 5, after a short transient in which a very steep 
drop occurs, the correlations decays much more slowly, 
indicating the long persistence of an average geometrical 
shape of the cluster. In the lower panel of Fig. [5] we com- 
pare the diffusion coefficient D(t) = (-^ J^iLi [ r i( r +0) ~ 
ri(/3)] 2 } of the 4 He atoms, as calculated from the original 
configurations of each quantum path generated during 
the simulation, with that obtained after rotating back- 
ward the configurations at time r in order to minimize 
the particles diffusion with respect to the configuration 
at time (3. The high degree of rigidity of the cluster struc- 
ture near the impurity is confirmed by the considerable 
effectiveness of the diffusion minimization procedure. 

In order to visualize the intrinsic shape of the 4 He clus- 
ter, we have calculated the helium density distribution 
p(r) = (Xa=i ^( r — r 0) around the ion, by applying 
the mentioned backward-rotation procedure. The con- 
tribution to p(r) from each quantum path sampled by 
the PIGS simulation is then further rotated to minimize 
the diffusion of the particles centroids of different paths. 
Fig [3] shows the density distribution for Li + @Hc 8 . The 
solvent accumulates in eight regions roughly arranged in 
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FIG. 4: Radial density distribution function of the He atoms 
in Li + @He 70 (solid black line), Na + @He 70 (dashed line), 
Be+@He 70 (red), Mg+@He 70 (green), Ca+@He 70 (blue). 



two parallel squares, rotated by 7r/4 with respect to each 
other. Indeed, the lowest non-vanishing values for a 
polyhedron with such a structure correspond to L = 4 
and L = 5. The results reported in Figs. [2] and [3] seem 
to indicate that Li + @Hes behaves much like a rigid body 
rather than as an ion embedded within a liquid shell. 



B. Multiple-shell ion-doped 4 He clusters 

We come now to the main concern of this paper, i.e. 
to discriminate between solid-like versus liquid-like struc- 
tures around helium-solvated ions. To this aim we have 
tried to apply the criterion of Ref. [HI to the more gen- 
eral case of clusters consisting of more than one solvation 
shell. 

Instead of extending the integral in Eq.rj]to the entire 
cluster, one may calculate multipole moments Qlm of a 
region within a suitably chosen distance from the center 
of mass. Thus, the criterion described in the previous 
Section allows to identify a solid-like behavior of that 
sub-region of the cluster. In Fig. [5] we show the radial 
solvent density profiles for 70 4 He atoms clusters doped 
with different ions. By taking the first minima of the 
curves as cut-off distances in the integral which defines 
Qlm, we accessed the imaginary-time auto-correlation 
functions of the multipoles of the first shell of the cor- 
responding clusters. Table I reports the integral of the 
4 He density within the first solvation shell for the various 
ions that we have considered. 



TABLE I: Integral of the 4 He density distribution within the 
first solvation shell of 4 He 7 o clusters doped with different ions. 



Li+ 


Na+ 


Be+ 


Mg+ 


Ca+ 


8.24 


12.02 


14.55 


19.02 


23.18 
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FIG. 5: Comparison of the multipole time-correlation func- 
tions cl (t) for L = 4 (upper panel) and L = 5 (lower panel) , 
as calculated for the whole Li + @He7o (empty circles), within 
the first solvent shell of the same cluster (filled circles), and 
for CO@Hei 5 (triangles). 

1. Alkali ions doped 4 He clusters 

In Fig. [5] we compare the multipole time-correlation 
functions cl(t) for L = 4 (upper panel) and L = 5 
(lower panel), as calculated for the first solvation shell in 
Li + @He7o, for the same whole cluster, and for CO@Hei5. 
The time correlations of the L = 4 and L = 5 multi- 
poles of the first shell, albeit differing in magnitude, both 
show a very slow decay (almost a saturation in the case 
of L = 4), whereas the same quantities, when considered 
for the whole cluster, decay rapidly to zero, very similarly 
to what happens in CO@Hei5. Also the L — 6 correla- 
tion function (not reported) exhibits a slow decay, just 
slightly more rapid than for the L = 4 and L = 5 curves. 
On the contrary, the decay of all the other cl(t) correla- 
tions, calculated within the first shell, is extremely fast. 
We interpret this result as an indication that the 4 He 
structure in the first solvation shell maintains an average 
shape, behaving like a solid cage which separates the Li + 
ion from the liquid environment of the outer shells. The 
signature of the rigidity of the first shell is lost among 
the random fluctuations in the more external regions, 
and thus it disappears in the multipole correlations of 
the entire system. 

At least qualitatively, we can attempt a better under- 



standing of the behavior of the time-correlations cl(t). 
In a very simplified and ideal picture, the diffusion of the 
4 He atoms of a solid first solvation shell during the sim- 
ulation can be thought as the combination of two kind of 
motions: one approximatively corresponding to a coher- 
ent rotation of the whole shell, and the other due to high- 
frequency small-amplitude fluctuations about the smooth 
trajectory of the former. When the cluster configurations 
along the path are rotated backward in order to mini- 
mize the particles diffusion, the neat rotation is totally 
suppressed, whereas the fast fluctuations are not. Thus 
the only surviving motions can be seen as the incoherent 
fluctuations of 4 He atoms in the first shell about average 
positions. More realistically, other motions, not removed 
by the backward rotation procedure, enter the decompo- 
sition of the 4 He atoms diffusion in the first shell. If these 
motions are sufficiently slow, i.e. if their energies are suf- 
ficiently low, the average shape of the first shell is not 
disrupted too quickly. Correspondingly, since imaginary- 
time correlation functions can be written as a sum of 
decaying exponentials whose decaying constants are the 
excitation energies of the system, the cl(t) correlations, 
characteristic of the system geometry, exhibit a slow de- 
cay 

Within the same qualitative picture, the atoms in the 
outer regions undergo random motions, none of which is 
reducible to a neat rotation. The slowest modes of the 
whole cluster are thus fast enough to prevent the persis- 
tence of any global shape (and to produce a rapid decay of 
the corresponding time-correlations cl(t)), but anyway 
they are much slower than the vibrations of the atoms in 
the first shell, thus explaining the lower decaying rate of 
the cl(t) correlation functions at small t with respect to 
those of the first shell. 

Fig. [5] depicts the 4 He density distribution within the 
radius of the first shell, where we found an integrated 
density of ~ 8.24 4 He atoms. The solvent concentrates 
mainly in eight regions arranged in two parallel squares 
rotated by tt/4 with respect to each other. At the oppo- 
site poles of this structure, along an axis perpendicular 
to the planes of the squares, two small lower-density ac- 
cumulations, probably due to incursions of He density 
from the second shell, are found (they disappear well be- 
fore the other eight ones when we increase the isodensity 
value chosen to draw the picture of Fig. [5]). As we ex- 
plained in discussing Li + @Hes, for a two rotated-squares 
structure the lowest non-vanishing qr, values correspond 
to L = 4 and L = 5. We note that the addiction of two 
more vertexes on a line perpendicular to the squares pro- 
duces a polyhedron having the lowest non- vanishing 's 
for L = 5 and L = 6. The behavior of the cl(t) correla- 
tion functions can thus be seen as the result of the con- 
comitant signals corresponding to these two structures. 
We stress at this point the need of considering dynamical 
correlations cl(t) of the multipoles of the cluster mass 
density, rather than their static magnitudes qr,, in order 
to reveal a solid-like behavior of the solvent. In fact, by 
calculating the quantities qL of the He density of the first 
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FIG. 6: He density distribution within the first solvation shell 
of Li + @He7o. The displayed isosurface contains 60% of the 
number of He atoms in the first shell. The Li + ion is located 
at the center of the picture, but it is not displayed. 

solvation shell of Li + @Hc7o, we obtain value at L = 4 
which is larger but of the same order of magnitude of 
the value at L = 3, whose corresponding imaginary-time 
correlation function, instead, decays extremely rapidly. 

Even more evident is the existence of a solid first shell 
in clusters doped with a Na cation. Our simulations indi- 
cate that the time correlations of the first five multipolcs 
within the first shell decay rapidly, whereas the L = 6 
shows an extremely slow decay. These results are com- 
patible with an icosahedral structure, whose lowest non- 
vanishing multipolcs correspond to L = 6 and L = 10. 
For the first shell of Na + @He7o, we show in Fig. [7] the 
slow-decaying L — 6 multipole time-correlation function 
(upper panel) and the 4 He density distribution (lower 
panel). The latter clearly shows an icosahedral order, in 
agreement with previous QMC calculations, based on the 
shadow wave functions technique iii 

Interestingly enough, despite the stronger interaction 
potential characterizing the Li + -Hc pair, we find a more 
pronounced solid-like behavior of the first shell of 4 He 
around a Na + ion than in the case of Li + . This is due to 
a subtle balance between the depth and the position of 
the minimum of the He-ion interaction potential, which 
determines the average number of 4 He atoms in the first 
solvation shell, and, hence, the average distances among 
them. For instance, the nearest neighbor He atoms in the 
first shell of Li + @He7o are separated by roughly 2.6 A, i.e. 
they find each other in the repulsive region of the He-He 
interaction potential. In contrast, the nearest neighbor 
distance in the first shell of Na + @He7o is slightly less 
than 2.8 A, that is within the well of the He-He pair 
potential and closer to its minimum (~ 3.0 A). 

The results for Na + @He7o have been obtained using 
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FIG. 7: Upper panel: multipole time-correlation functions 
cl(t) for L — 6, as calculated for the whole Na + @He7o clus- 
ter (empty circles) and within the first solvent shell (filled 
circles). Lower panel: He density distribution in the first 
shell of Na + @He7o. The isodensity surface contains the 60% 
of the density of the first shell. The Na + ion resides at the 
center of the structure and it is not displayed. 



a time step e = 0.0125 K _1 . By integrating the 4 He 
density within the first shell of Na + @He7o, we obtain an 
estimate of twelve 4 He atoms, in agreement with shadow 
wave function and path-integral Monte Carlo (PIMC) 
results ; 11 ! 27 but at variance with other existing PIMC 
calculations fi 2 - Our result is stable, within the statisti- 
cal error, when reducing the time step by a factor 2.5, 
thus showing a negligible time-step dependence in this 
e range. For larger time steps, however, the reliability 
of the simulations is not guaranteed. In Ref. [IH, where 
a time step of 1/20 K _1 is used, sixteen 4 He atoms are 
found, instead of twelve, in the first shell of a Na + @Heioo 
cluster at the temperature of 1 K. The above discrepancy 
may originate from the choice of the Hc-Na + interaction 
potentiate or, more likely, it is a time-step effect. In 
fact, for strongly attractive potentials, such as the Li + - 
He and the Na + -He ones, we noted that the use of a 
large time step (like the one used in Reffl^l may lead to 
overestimate the 4 He density in the first shell. In par- 
ticular, by simulating Na + @He7o with the same rather- 
large time step e = 0.05K -1 , we reproduced the result 
of Ref. [l^. We also observed that the 4 He density in 
the first solvation shell is much more dclocalizcd, lacking 
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FIG. 8: Multipole time-correlation functions cl(t) for L = 7, 
as calculated for the whole Be + @He7o cluster (empty circles) 
and within the first solvation shell (filled circles). 



the icosahedral structure previously found, and that the 
lowest twelve multipolcs time-correlation functions decay 
now much more rapidly. 



FIG. 9: He density distribution in the first shell of Be + @He7o. 
The isodensity surface contains the 60% of the density of the 
first shell. The Be + ion resides at the center of the structure 
and it is not displayed. 



2. Alkaline- earth ions doped 4 He clusters 

We now turn our attention to the case of 4 He clus- 
ters doped with alkaline-earth ions. Due to their larger 
radii and their relatively weaker interactions with helium 
atoms, these ions are expected to give rise to bubble-like 
structures of the 4 He liquid around the impurity. This 
expectation is confirmed by mobility measurements at fi- 
nite temperature for several alkaline-earth ions but not 
for Bc + , for which results of similar measurements sug- 
gest instead a snowball structure of the defect.— For this 
reason, in the following, we will focus particularly on the 
solvent structure around Bc + . The rigidity of the defect 
structures of Be + and Mg + can hardly be distinguished 
on the basis of static correlation functions among 4 He 
atoms^ 4 - thus spurring us to consider the multipole mo- 
ments dynamical correlations criterion.— 

In order to detect even weak signatures of a possible 
solid-like order around the alkaline-earth ions, the clus- 
ters configurations, used for calculating both static and 
dynamical quantities, have been rotated back in order to 
minimize the particles diffusion within the first solvation 
shell, and not in the whole clusters. According to our 
results, all the multipolcs correlations of the first shell of 
Be + @He7o up to L = 5 decay rapidly (even faster than 
in the whole cluster), whereas for L = 6 and, more ev- 
idently, for L = 7 the correlation decay is much slower 
than that of the whole cluster. Also in this case the mag- 
nitude of the multipoles q^ provide a guide in the search 
of the slow decaying time-correlations. Fig. [5] displays 
the L = 7 multipole correlation in the first solvation shell 
and in the whole cluster. In spite of the faster decay than 
for Li + @He 70 and Na + @Hc 70 , the behavior of Cj,(t), for 
L = 7, indicates an extremely slow disappearance of the 



shape assumed by the solvent around the Be + ion. The 
average shape of the solvent in the first shell is shown 
in Fig. [9] Despite being less localized than for alkali ion 
doped clusters, the 4 He density in the vicinity of Bc + 
appears anyway significantly structured. 

Fig. [TU] compares the diffusion coefficient in imaginary 
time, D(t), as calculated after the backward rotation, 
for the 4 He atom in the whole cluster and within the 
first shell, for Li + @He 70 and Be+@He7 . For both sys- 
tems, with increasing the (imaginary) time r, the coeffi- 
cient D(t) of the whole cluster increases much more than 
within the first solvation shell, due to the liquid character 
of the external regions. We interpret the small diffusion 
rate of the He atoms in the first shell, and, moreover, the 
close similarities in D(t) for Be + @He 7 o and Li + @Hc 7 o 
(which we consider a prototype of a cluster with a first 
shell with solid-like order) as an indication of a nearly 
rigid-body behavior of the first shell of Be + @He 7 o. This 
conclusion will be strengthened by the comparison with 
the results for Mg + @He7o, which will be discussed next. 

For Mg + @He7o the calculation of qL, within the first 
solvation shell, gives a faint signal for L = 8. How- 
ever, the corresponding imaginary-time correlation func- 
tion cl(t), which is reported in Fig. QT1 exhibits an ex- 
tremely fast decay, despite being the slowest decaying 
one among those we have analyzed (up to L = 12). Our 
findings suggest that the He atoms surrounding Mg + un- 
dergo random fluctuations which disrupt the local struc- 
ture of the solvent much more quickly than in the other 
cases examined so far, thus preventing the formation of 
any longly persisting solid-like order. For Mg + @He 7 o the 
rapid decaying of the cl (t) correlations of the first shell 
is very similar to that of the corresponding ones of the 
entire cluster, thus indicating a liquid-like behavior of the 
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FIG. 10: Upper panel: diffusion coefficient D(r) as a function 
of (imaginary) time r calculated from the backward-rotated 
configurations, for the He atoms in the whole Li + @He7o clus- 
ter (empty circles), and for the He atoms in its first solva- 
tion shell (filled circles). Lower panel: same quantities for 
Be + @He7o, displayed with the same symbols. 

whole system. 

As it is shown in the upper panel of Fig. [T^l the 4 He 
density distribution in the first shell is considerably more 
diffuse than what found for the other ions. By increas- 
ing the isodensity value chosen for the graph, it becomes 
possible to identify regions of higher density, nevertheless 
the solvent is much more delocalized than for Be + @He7Q. 



1 f- 




FIG. 11: Multipole time-correlation function cl(t) for L = 8, 
as calculated within the first solvation shell of Mg + @He7o- 
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FIG. 12: Upper panel: He density distribution in the first 
shell of Mg + @He7o . The integral of the density contained 
within the displayed isodensity surface amounts to 60% of the 
number of He atoms in the first shell. Lower panel: diffusion 
coefficient D(t) as a function of (imaginary) time r for the 
He atoms in Mg + @He7o, as calculated form the backward- 
rotated configurations, for the whole cluster (empty circles) 
and for the first solvation shell (filled circles). 



The lower panel of Fig. [T2l compares the imaginary-time 
diffusion of the He atoms in the entire cluster and in the 
first solvation shell, after the backward rotation. The two 
curves are now closer and they even merge within the er- 
ror bars, thus proving that the rotational diffusion of the 
particles in the first shell cannot be reduced to a neat 
rotation. The liquid-like behavior of the solvent is even 
more evident in Ca + @He7o, for which the first shell mul- 
tipole time-correlation functions decay very rapidly, the 
density distribution in the first solvation shell is largely 
delocalized and the diffusion coefficients of 4 He atoms of 
the inner and of the outer shells are much less distin- 
guishable. 



IV. CONCLUSIONS 

Using PIGS simulations, we have studied the solvation 
of alkali and alkaline-earth ions in 4 He clusters. In or- 
der to discriminate between solid-like versus liquid-like 
behavior of the solvent in the vicinity of the impurity, 
we employed a criterion based on dynamical correlations 
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of the multipole moments of the 4 He density in the first 
shell, previously proposed for quantum clusters consist- 
ing of just one completed shell. We also used other indi- 
cators to clarify the structure of 4 He around the impurity. 

The dynamical correlations criterion, whose applicabil- 
ity to clusters with more than one shell was not granted 
a priori, proved suitable to address situations in which 
regions with different rigidity coexist within the system, 
wit hout requiring any prior knowledge of the defect struc- 
ture. In particular we find bubble (liquid-like) structures 
for heavier alkaline-earth ions (Mg + and Ca + ), whereas 
alkali-ions are caged by a solid-like, snowball structure, 
further embedded in an outer liquid environment. For 



these two categories, whose features are expected to be 
different both theoretically^! 4 - and experimentally^ the 
multipole correlations signatures are clearly different and 
unambiguously distinguishable. 

For Be + , our simulations indicate a border-line behav- 
ior, which, however, has much more similarities with al- 
kali ions than with heavier alkaline-earth ones. This re- 
sults suggests that Be + forms a snowball upon solvation, 
possibly elucidating the "anomalous" mobility of Be 4 " 
in liquid helium with respect to Mg + , Ca + and other 
alkaline-earth ions. 6 
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